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ABSTRACT 



Oh 



Context. Protostellar jets from intermediate- and high-mass protostars provide an excellent opportunity to understand the mechanisms 
responsible for intermediate- and high-mass star formation. A crucial question is if they are scaled-up versions of their low-mass 
counterparts. Such high-mass jets are relatively rare and, usually, they are distant and highly embedded in their parental clouds. The 
IRAS 20126+4104 molecular jet, driven by a 10 4 Lq protostar, represents a suitable target to investigate. 
Aims. We present here an extensive analysis of this protostellar jet, deriving the kinematical, dynamical, and physical conditions of 
j_( . the H2 gas along the flow. 

Methods. The jet has been investigated by means of near-IR H2 and [Fe 11] narrow-band imaging, high resolution spectroscopy of the 
1-0 S(l) line (2.12/mi), NIR (0.9-2.5 ^m) low resolution spectroscopy, along with ISO-SWS and LWS spectra (from 2.4 to 200 yum). 
Results. The flow shows a complex morphology. In addition to the large-scale jet precession presented in previous studies, we detect 
a small-scale wiggling close to the source, that may indicate the presence of a multiple system. The peak radial velocities of the H 2 
knots range from -42 to - 14 km s" 1 in the blue lobe, and from -8 to 47 km s -1 in the red lobe. The low resolution spectra are rich in H 2 
emission, and relatively faint [Fen] (NIR), [Oi] and [Cn] (FIR) emission is observed in the region close to the source. A warm H2 
gas component has an average excitation temperature that ranges between 2000 K and 2500 K. Additionally, the ISO-SWS spectrum 
reveals the presence of a cold component (520 K), that strongly contributes to the radiative cooling of the flow and plays a major 
role in the dynamics of the flow. The estimated of the jet is 8.2+0.7 Lq, suggesting that IRAS 20126+4104 has an accretion 
rate significantly increased compared to low-mass YSOs. This is also supported by the derived mass flux rate from the H2 lines 
(M,„,(H2)~7.5xlO~ 4 Mq yr~'). The comparison between the H? and the outflow parameters strongly indicates that the jet is driving, 
at least partially, the outflow. As already found for low-mass protostellar jets, the measured H 2 outflow luminosity is tightly related to 
the source bolometric luminosity. 

Conclusions. As for few other intermediate- and high-mass protostellar jets in the literature, we conclude that IRAS 20126+4104 jet 
is a scaled-up version of low-mass protostellar counterparts. 
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1 . Introduction mechanism of these latter objects is still debated. The observa- 
tions are, however, strongly limited by the large distance of the 
Protostellar jets and outflows are a ubiquitous phenomenon massive star forming regions, the considerable extinction, and 
among young stellar objects (YSOs) of different mass and lu- the short uf etime of maS sive YSOs. In addition, these objects 
minosity (see e. g. Shepherd 2003). They are usually explained are often grouped in smaU clusters, which confuses the mor- 
as a consequence of accret ion fro m a disc ar ound the protostar phology of massive star forming regions even more. Therefore, 
(seee.g.Pudntz&Norman[198aCamenzind[199g.Thisispar- optical and IR studies of intermediate- high-mass protostel- 
ticularly true for low- and partially, for intermediate- and high- lar jets are very rare ^ and only a few examples are present 
mass YSOs up to L bol ~10 4 L Q (or spectral type B0), where col- in literat ure: e.g. IRAS 18162-2048 {L bol ~2xlO 4 L , Martlet 
hmated outflows, often d riven by protostellar jets, have been ob- al [1995] !, tr A S 20126+4104 (L bol ~10 4 L q, Aya la et al. [1998 ), 
served (Shepherd 2003). On the contrary, no highly collimated IRAS 15547.4247 (6xl0 4 L o , Brooks et al.|2003), IRAS 18151- 
outflow or circumstellar disc have been observed in high-mass 1208 (2xlO 4 L , Da vis et al. 120041 and IRAS 1 1 101-5829 
protostars exceeding 10 5 Lq (or O-type stars, the spectral type (1Q 4 Lq; Gredel 2006k the M1? disc si]houette {Myso „ 15 Mq> 
and L boI o f theM 17 disc silhouette a re not clear yet) (see e. g. Narn berger et al. 120071 Expanding the number of observations 
Arce et al. 120071 Zinnecker & Yorke|2007J, and the formation of intermediate . and high-mass jets and comparing their gen- 

eral properties with those of low-mass protostellar jets is there- 

Send offprint requests to: A. Caratti o Garatti e-mail: fore important to understand if differences exist, or if high-mass 

caratti@tls-tautenburg.de protostellar jets are just scaled up versions of their low-mass 

* Based on observations collected at UKIRT, TNG, and at Subaru and counterparts. We therefore investigate the kinematical and phys- 
ISO data archive. 
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ical properties of the IRAS 20126+4104 jet by means of NIR 
narrow-band imaging, high resolution and low resolution IR 
spectroscopy. We then compare our findings with those of other 
high- and low-mass protostellar jets available in literature. 

IRAS 20126+4104, at a distance 1.7kpc, i s a ver y-well stud- 
ied high-mass YSO (M ~7 Mq, Cesaroni et al. |1997l Cesaroni et 
al. 120051 ), in a very early stage of evolution. It is accreting mass 
at a very high rate (M acc ~2xl0~ 3 M©yr _1 , Cesaroni et al. 2005 ) 
and it gives birth to a large poorly-collimated CO outflow. It har- 
bours the first H2 jet detected from a high mass YSO (Ayala et 
al. II 9981 ), which had previously been seen in SiO emission close 
to the source (Cesaroni et al. 1999 ). The H2 jet extends for about 
1 pc. Its 'S' shape morphology suggests that it is precessing with 
a period of ~60 000 yr and with a wide precession angle of about 
37° (Shepherd et al. 2000]). As a consequence, the inclination of 
the flow with respect to the plane of the sky changes strongly 
from ~9°, close to the source, up to ~45° in the outer part of the 
flow. 

The structure of this paper is as follows. In Sect.|2]our obser- 
vations are presented. Sect. [3]reports an overview of our results, 
including the morphology of the H2 jet, its kinematics, the physi- 
cal parameters of the gas, its energy and mass flux rate. In Sect. [4] 
we briefly consider the cause of the newly detected small-scale 
precession mode. Then, we discuss our H2 data in relation to 
the CO outflow literature data. Finally, we compare the proper- 
ties of the IRAS 20126+4104 jet with other high- and low-mass 
protostellar jets. 



2. Observations and data reduction 

Our data were collected at the UK Infrared Telescope (UKIRT), 
and at the 3.5-m Italian Telescopio Nazionale Galileo (TNG). 
Further data have been retrieved from the Subaru and ISO 
archive^ The relevant information on the observational settings 
is summarised in Table Q] 



2.1. Imaging: H 2 &[Fe\\] 

We used narrow-band filters centred on the H2 (2.12/im) and 
[Fen] (1.64yum) lines to detect both molecular and ionic emis- 
sion along the flow. Additional broad-band images (H and K') 
were gathered to remove the continuum. Our images were col- 
lected at UKIRT, using the near-IR instrument UIST (Ramsay 
Howat et a l. 12004] ) ([Fen], H), and at the TNG, using NICS 
(Baffa et al. 200TJ (H2, K'). All the raw data were reduced using 
IRAF0 packages, applying standard procedures for sky subtrac- 
tion, dome flat-fielding, bad pixel and cosmic ray removal, and 
image mosaicing. The resulting mosaics cover a region around 
IRAS 20126+4104 of 3:4x3:4 for UIST, and 5:6x5:6 for NICS. 
The calibration for both instruments was obtained by means of 
photometric standard stars observed in both narrow- and broad- 
band filters. In the calibrated and continuum-subtracted narrow- 
band images, we measured H2 and [Fen] fluxes for each de- 
tected knot using the task polyphot in IRAF, defining each region 
within a 3<x contour level above the sky background. 

Additional high-resolution narrow-band images (H2, Br r , 
and Kcont) of the central part of the molecular flow (~35"x35") 



1 These data are available at http://smoka.nao.ac.jp/ (Subaru), and 
http://www.iso.vilspa.esa.es/ida/index.html (ISO) 

2 IRAF (Image Reduction and Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are operated 
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Fig. 1. Close up of the TNG H2 mosaic, showing the entire 
flow with superimposed positions of the slits, and the ISO- 
SWS FoV. The positions of the H2 knots and the IRAS source 
(a = 20' ! 14 m 26.03 ! and 6 =41°13'32".58 [12000], Hofner et 
al. 2007) are also indicated. 
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were retrieved from the Subaru data archive (Baba et al. 2002). 
The set of data were taken with the Coronagraphic Imager with 
Adaptive Optics (CIAO) (Murakawa et al. 2003 ), used as a high- 
resolution NIR imager (07022/pixel). These data were used for 
a morphological study of the flow close to the source only (see 

Sect.OandgZQ. 

Finally, we used K-band stars from the Two Micron All Sky 
Survey (2MASS) for the astrometric calibration of all the im- 
ages. As a result the formal errors for the plate solutions are 
-075, 074, and 072, for the NICS, UIST, and CIAO images, re- 
spectively. 

2.2. Spectroscopy 

2.2.1. Low resolution spectroscopy 

NIR low resolution spectroscopy (H ~500, slit width 1") was 
acquired with NICS during two runs at the TNG, using three 
different grisms covering all the full NIR spectrum (0.9-2.5/im) 
(see Table [T]). The position of the slits is shown in Fig.Q] Slits 
1 and 3 are positioned on the knots close to the source (knots 
A, B, and C) and parallel to the jet axis, covering almost all the 
H2 emitting area. Slits 2 and 4 encompass the so called 'jet' (see 
Ayala et al. 1 19981 ), that is the farthest H2 emission in the red lobe 
(see also Fig.Q~|i. 

To perform our spectroscopic measurements, we adopted the 
usual ABB'A' configuration, with a total integration time be- 
tween 1800 s and 2400 s per grism per slit. Each observation 
was flat fielded, sky subtracted and corrected for the curva- 
ture caused by longslit spectroscopy, while atmospheric features 
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IMAGING 



Date of obs. 


Telescope/ 


Filter 


Resolution 


seeing 


Exp. Time 


(d,m,y) 


Instrument 


Band 


("/pixel) 


(") 


(s) 


01.10.2006 


UKTRT/UIST 


[Fen] 


0.12 


0.6 


720 


01.10.2006 


UKTRT/UIST 


H 


0.12 


0.6 


360 


06.08.2006 


TNG/NICS 


H 2 


0.25 


1.2 


900 


06.08.2006 


TNG/NICS 


K' 


0.25 


1.2 


90 


10.07.2003 


Subaru/CIAO 


H 2 , Bry, K t . „, 


0.022 


0.9 


1800 



SPECTROSCOPY 



Date of obs 


Telescope/ 


. Wavelength 


Hiti 


PA. 


Slit/Aperture width 


K 


Encompassed Knots 


Notes 


(d,m,y) 


Instrument 


(fim) 


(s) 


C) 


(") 








01.10.2006 


UKIRT/CGS4 


2.1218 


2700 


300 


0.5 


18 500 


A1,A2,X,B,C2 


UKIRT-Slitl 


01.10.2006 


UKIRT/CGS4 


2.1218 


2700 


300 


0.5 


18 500 


A1,A2,B,C1 


UKIRT-Slit2 


01.10.2006 


UKIRT/CGS4 


2.1218 


2700 


347 


0.5 


18 500 


jet - knots 1,2,3,4 


UKIRT-Slit3 


01.10.2006 


UKIRT/CGS4 


2.1218 


2700 


347 


0.5 


18 500 


D,C2 


UKIRT-Slit4 


06.08.2006 


TNG/NICS 


0.88-1.45 


2000 


299 




500 


A,B,C 


TNG-Slitl 


06.08.2006 


TNG/NICS 


1.09-1.80 


1800 


299 




500 


A,B,C 


TNG-Slitl 


06.08.2006 


TNG/NICS 


1.40-2.47 


1800 


299 




500 


A,B,C 


TNG-Slitl 


06.08.2006 


TNG/NICS 


1.40-2.47 


1800 


350 




500 


jet 


TNG-Slit2 


18.10.2003 


TNG/NICS 


1.09-1.80 


2400 


300 




500 


A,B,C 


TNG-Slit3 


18.10.2003 


TNG/NICS 


1.40-2.47 


2400 


300 




500 


A,B,C 


TNG-Slit3 


18.10.2003 


TNG/NICS 


1.09-1.80 


1800 


353.5 




500 


jet 


TNG-Slit4 


18.10.2003 


TNG/NICS 


1.40-2.47 


2400 


353.5 




500 


jet 


TNG-Slit4 


11.05.1996 


ISO/SWS 


2.38^15.2 


1912 




33x20 


1000-2000 


A,B,C,D 




30.12.1995 


ISO/LWS 


43-197 


750 




80 


200 


A,B,C,D 





were removed by dividing each spectrum by a telluric standard 
star (spectral O type). The wavelength calibration was obtained 
from Xenon and Argon lamps. The flux was calibrated observing 
the photometric standard star AS 35-0 (Hunt et al. 119981 ). As the 
used grisms overlap between ~1.1 and 1.8 fan (see Tab. [TJ, the 
single spectra were combined to obtain deeper spectra. Finally, 
since no substantial difference was found between the pairs of 
spectra from Slit 1 and 3, nor from Slit 2 and 4, they were com- 
bined together to get a deeper spectrum of each observed knot. 

2.2.2. High resolution spectroscopy 

Our H 2 1-0S(1) (A mc = 2.1218356 /mi; Bragg et al. [19821 
echelle spectra were obtained at UKIRT using the spectrome- 
ter CGS4 (Mountain et al. 1 19901 1 (see Table Q}, equipped with a 
256x256 pixel InSb array (0741 xO^O pixel scale, in the disper- 
sion and spatial direction, respectively). A 1 -pixel-wide slit was 
used, resulting in a velocity resolution of ~7.4k.ms~ 1 pixel -1 , 
and corresponding to a nominal resolving power of ~ 1 8 500. The 
instrumental profile in the dispersion direction, measured from 
Gaussian fits to the OH sky lines, was -10.9 km s _1 (or 1.46 A). 

The targets and position angles (P.A.s) of the four slits are 
reported in Table [2 while the slit positions are shown in Fig.Q] 
and encompass all the H2 knots of the flow. Slit 1 and 2 are 
parallel to the flow axis and positioned on the knots close to the 
source (knots A, B, and C). Slit 3 encompasses the so-called 
'jet', and Slit 4 encompasses knot D and partially knot C. The 
total exposure time for each position is 2700 s. 

The raw spectra were wavelength calibrated using the bright 
OH lines (Rousselot et al. 2000) detected on each frame. The 
employed IRAF tasks (identify, re-identify, fitcoords, and trans- 
form) also correct for the spatial distortion, and a 3rd-order fit 



in two dimensions has been used. As a result, the calibration is 
accurate to ~3kms _1 . The spectra were not flux calibrated. 

2.2.3. ISO Data Archive 

Mid-IR and far-IR spectra of the flow were retrieved from the 
ISO satellite Data Archive, in the form of Highly Processed Data 
Product (HPDP set 35500738, by W. Frieswijk et al.; HPDP set 
04300333 by Lloyd C, Lerate M., and Grundy T). The first set 
of data (from 2.38 to 45.2 yum, see also Tab.Q} was taken with the 
Short Wavelength Spectrometer (SWS, de Graauw et al. |1996b in 
the AOT01 grating mode ("R -1000-2000), and is centred on the 
IRAS source covering a region of 33"x20". The covered FoV is 
shown in Fig.Q] The second set of data (from 43 to 197 fim) was 
taken with the Long Wavelength Spectrometer (LWS, Clegg et 
al. [19961 ) in the AOT01 grating mode (A/AA ~ 200), and is po- 
sitioned ~48" north of the source. With a beam of 80" it covers 
large parts of the flow, except the so called 'jet'. 



3. Results 

3.1. H 2 and [Fen] imaging 

An overview of the large-scale morphology of the H2 flow is 
given in Fig.Q] The emission in the south-east, the so called 'jet', 
appears to be composed of four different knots, labelled knots 1 
to 4 (see Sect. 13.3b . No further knots were detected SE or NW 
up to ~2'.S (or ~1.4pc) from the source, down to a 3<x limit of 
~8x 10~ 17 erg s _1 crrT 2 arcsec~ 2 . 

However a new emission feature, labelled X 
(or(72000)=20'T4 ffl 26.1 ! and 5(72000)=41°13'32".4), has 
been detected close to the position of the IRAS source. Fig. [2] 
shows a close-up in H2 of the source position and the new knot, 
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2 1 0-1-2-3 -4 
R.A. Offset (arcsec) 



Fig. 2. Close-up of the H2 image around the presumed source 
position (Hofner et al. 120071 see also the caption of Fig.[TJ indi- 
cated by a circle (the radius represents the formal error for the 
NICS plate solution). [Fe 11] contours (5 in red, 30 in yellow, 50, 
60, 70 ,80 ,100 xcr ~ 2.5x 10~ 17 erg s -1 cirT 2 arcsec -2 ) are over- 
laid on the image. The positions of the new knot X, knot B, the 
source, and two field stars are also indicated. 



located in the redshifted lobe. The circle represents the formal 
error for the NICS plate solution centred on the IRAS source 
position (Hofner et al. 2007). Contours of the [Fen] emission 
are superimposed on this image as well. The two stars in the 
field have been used to match the H2 and [Fe n] images exactly. 
Both our H2-K' and [Fen]-H data (as well as the spectroscopic 
data, see Sect 13. 2\ indicate that the emission of this knot arises 
from a shock. Remarkably, the knot is overlapping with one 
of the bipolar features (located SE of the source) detected 
in the K-band by Sridharan et al. (2005). No further [Fen] 
emission has been detected along this flow, down to a 3<x limit 
of ~3x 10 -16 erg s _1 cirT 2 arcsec -2 , nor along the blue lobe at 
the position where a Ha and [S 11] emitting knot was observed 
by Shepherd et al. j2000] l. 

The higher angular resolution of Subaru images delineates 
details of the morphology of the H2 flow close to the source. In 
Fig. [3] we show the continuum subtracted H2 image. The Bry 
image turned out more efficient than the K co „, image in remov- 
ing the strong continuum detected on the knots. Its use is justi- 
fied by the absence of any Bry emission in the observed knots 
(see Sect. 13.21 1. Scattered H2 emission would, of course, not be 
subtracted by this procedure. The H2 flow is wiggling notice- 
ably with an opening angle of -16° around an axis with RA. of 
— 61°. This fully agrees with SiO radial velocity observations 
(Cesaroni et al. 119991 Su et al. |2007l l, that revealed a precessing 
jet with an opening cone of 21°, RA — 60°, and an inclination of 
9° with respect to the plane of the sky. This small-amplitude pre- 



41:13:40.0 




37.8 




35.6 











33.5 


<n 




^5 




u 


31.3 






a 






29.2 




27.0 



H-Br Y 



star - 9| 
star - • 

2 1 3 



star - « 

O 1 

Source 



B 



star 



27.0 



26.5 



20:14:26.0 
R.A. (J2000) 



25.5 



Fig. 3. Continuum- subtracted H2-Bry image from Subaru. The 
Bry image was used as continuum, since no Bry emission is de- 
tected on our spectra (see text). The identified knots and the pre- 
sumed source position (see also the captions of Fig.Q]and Fig.O 
are indicated. At this scale the wiggling structure of the flow be- 
comes visible. 



cession seems to be superimposed to the larger one observed at 
large scale (Shepherd et al. 2000 Cesaroni et al. 2005 see also 
Fig- ID- In the blue lobe (see Fig. [3]) knot B appears composed of 
a faint jet (feature 1) ending in a bow-shock (2), that is preceeded 
by two smaller bow-shocks (3). Further ahead, knot C consists 
of at least two bright structures. In the red lobe, the new knot X 
is visible as well, and, farther out, knot A appears fragmented in 
2-3 sub-structures. 

3.2. Low resolution spectroscopy 
3.2.1. Low resolution - TNG 

From our four TNG slits, we have derived low resolution spectra 
of six knots, namely A, B, C, and X close to the source, and knot 
1 and 4 along the 'jet' . All these spectra are rich in H2 emission. 
Remarkably, no ionic emission is detected along the flow, with 
the exception of faint [Fe 11] emission in knot X (also observed in 
the narrow -band image), and in knot B (below the detection limit 
of our imaging). It is worth to note that the ratio between the 
H2 lines (see, e.g., 2.12/im) and the [Fen] lines changes in the 
observed knots. Continuum emission from the regions enclosing 
knots A and B is observed in the K band. 

The lines detected in the spectra, along with their vacuum 
wavelengths, and fluxes (uncorrected for the extinction), are pre- 
sented in Table [2] Line fluxes have been obtained by fitting the 
profile with a single or double Gaussian in case of blending. The 
uncertainties associated with these data derive only from the rms 
of the local baseline multiplied by the linewidths. Lines showing 
fluxes with a S/N ratio between 2 and 3 have been labelled in the 
Table. 

Figure [4] shows the spectrum of knot C, the brightest one ob- 
served. Only H2 lines are visible from 1 to 2.5 fim, involving vi- 
brational levels from v=l up to v=5 (coming from energy levels 
from 6000 K up to 30 000 K). The presence of lines from high-v 
levels (i. e. with high excitation energies) indicates a high exci- 
tation of the gas. In the remaining knots we identify H2 emission 
lines up to v<3, mostly detected in the H- and K-bands (except 
for a few v=2 lines from the J-band observed in knot B). 

Considering both the combined effects of the knot morpholo- 
gies and the slit positions, there is an excellent agreement be- 
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Table 3. Observed lines in the IRAS 20126+4104 outflow: ISO 
lines 



Term ifrim) F ± A F(l(T' 3 erg cm' 2 s"' ) 

H2 Lines Flux 

1-0 Q(l) 2.407 1.8±0.4 

1-0 Q(2) 2.413 1.4±0.4 

1-0 Q(3) 2.424 0.8±0.4 

1-0 Q(4) 2.437 0.7±0.3: 

1-0 Q(5) 2.454 1.4±0.3 

1-0 Q(6) 2.476 1.3±0.2 

1-0 Q(7) 2.500 1.4±0.1 

1-0 0(2) 2.627 0.6±0.1 

1-0 0(3) 2.803 0.7±0.1 

1-0 0(4) 3.004 0.6±0.2 

1-0 0(5) 3.235 0.8±0.1 

0-0 S(6) 6.109 21±3 

0-0 S(5) 6.909 21 ±3 

0-0 S(4) 8.026 30±3 

0-0 S(3) 9.665 32±6 

0-0 S(2) 12.278 26±8 

0-0 S(l) 17.034 56±8 

ionic lines 

[Oi]a 3 Pi -c?P 2 63.18 260±10 

[Cn]a 3 P 3/ 2 -a 3 Pi/2 157.74 370±10 




wavelength {fJ.m) 



Fig. 5. ISO- SWS and LWS detected lines. 

tween the photometric fluxes in our images and the fluxes mea- 
sured in our spectra for both [Fen] (1.64 /mi) and H2 (2.12/mi) 
lines. 



Table 4. H2 radial velocities of individual knots in the 
IRAS 20126+4104 flow from CGS4 observations. 



CGS4 slit 1 



knot 


Vrad 


FWZI 




(kms- 1 ) 


(km s- 1 ) 


Al 


2±3 


110 


A2 


9±3 


150 


X 


-8±3; -28±3; 53±3 


175 


B 


-14±3 


160 


C2 


-25±3 


160 


CGS4 slit 2 


Al 


-1±3; -52±3 


130 


A2 


9±3;-12±3;53±3 


180 


B 


-14±3; -46±3; 26±3; 64±3 


170 


CI 


-19±3 


130 


CGS4 slit 3 


jet knot 1 


47±3; 19±3 


70 


jet knot 2 


27±3 


75 


jet knot 3 


41+3; 2±3 


90 


jet knot 4 


6±3; 36±3 


90 


CGS4 slit 4 


D 


-42±3; -17±3 


70 


C2 


-28±3; -5±3; 3±3 


90 



are detected, superimposed on a strong continuum observed in 
both SWS and LWS spectra. The detected lines are shown in 
Fig. H2 lines from the v=l and v=0 levels have been de- 
tected between 2.4 and 17 fim. In particular, the 1-0Q lines (also 
observed in our NIR spectra), the 1-0 O lines, and the 0-0 S 
lines, that usually trace the coldest component of the H2 flow 
(see e. g., Froebrich et al. 120021 Giannini et al. 2004 Giannini 
et al. 120061 1 have been detected. In addition, two bright ionic 
lines have been observed in the LWS spectrum, namely [0 1] at 
63 /im and [Cn] at 158//m. The origin of these lines is indeed 
intriguing, because no strong ionic emission lines were detected 
in the NIR. Possibly, they could originate from the same spot 
where [Fe 11] is detected, and/or in a deeply embedded region 
(close to the source) not visible at NIR wavelengths, because 
our spectra do not show any evidence of the bright [C 1] dou- 
blet at 1 yum, often detected in low-mass jets (see e. g. Nisini et 
al. 120021 Giannini et al. 2004). However, a more convincing ex- 
planation is that these lines are not arising from the shocks, but 
from a photo-dissociation region (PDR) around the massive ob- 
ject. Indeed, from the [Oi](63^m)/[Cn](158//m) ratio, we can 
discriminate between the two scenarios (see e. g., Hollenbach & 
McKee 1989). For intensity ratios below 10 the PDR origin is 
favoured. In our LWS spectrum this ratio is close to 1 . 



3.2.2. Low resolution - ISO 

Table |3]lists the lines identified in the ISO-SWS and ISO-LWS 
spectra, their vacuum wavelengths and integrated fluxes, esti- 
mated from Gaussian fits to the (unresolved) line profiles, to- 
gether with the error from the rms noise of the local baseline. 
Due to the large FoV of the two instruments, it is not possible 
to determine the spatial extension of the emitting region, thus 
the measured fluxes can originate from any or all of the dif- 
ferent knots (A- • -D) close to the source. Several emssion lines 



3.3. High resolution spectroscopy 

In Table [4] we report the results of the high resolution H2 spec- 
troscopy. For each knot in each slit we indicate the H2 radial ve- 
locities (vrad), corrected for the cloud speed with respect to the 
LSR {vlsr = -3.5kms _I , Cesaroni et al. 1997). Two or more 
velocity components are given where detected. The first reported 
value is the peak value. The last column gives the full width at 
zero intensity (FWZI) of the line profile, measured where the flux 
reaches a 2<x background noise level (see e. g. Davis et al. 1200 11 1. 
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Table 2. Observed lines in the IRAS 20126+4104 outflow: knots A, B, C, X, jet-1, jet-4. (To be continued) 



Term 


Aifim) 




F ± A F(10- 15 ergcrrr- s" 1 ) 


H2 Lines 


A 


B 


C X knot 1-jet knot 4-jet 


2-0 S(9) 


1.053 




1.3 t 0.1 


2-0 S(S) 


1.057 




0.9, tU 


2-0 S(7) 


1.064 




2.8±0.4 


2-0 S(6) 


1.073 




1.6±0.4 


2-0 S(5) 


1.085 




3.5±0.4 


2-0 S(4) 


1.100 




1.8+0.3 


2-0 S(3) 


1.117 


1.3±0.4 


4.3±0.2 


3-1 S(9)+3-l S(10) 


1.120 




2.0±0.2 


3-1 SC11) 


1.122 




1.0±0.3 


3-1 S(8) 


1.125 




0.7+0.2 


3-1 S(7) 


1.130 




1.4±0.2 


3-1 S(13) 


1.132 




0.5±0.2: 


2-0 S(2) 


1.138 




2.0±0.2 


3-1 S(6) 


1.140 




0.7±0.2 


3-1 S(5) 


1.152 




1.9+0.3 


2-0 S(l) 


1.162 


1.3±0.4 


3.5±0.3 


3-1 S(4) 


1.167 




1.2±0.3 


3-1 S(3) 


1.186 




2.5±0.3 


2-0 S(0)+4-2 S(10) 


1.189-1.190 




0.9±0.2 


4-2 S(9)+4-2 S(8) 


1.196-1.199 




0.7+0.2 


4-2 S(7) 


1.205 




1.1+0.2 


3-1 S(2) 


1.207 




0.9+0.2 


4-2 S(6) 


1.214 




0.6+0.2 


4-2 S(5) 


1.226 




1.4+0.2 


3-1 S(l) 


1.233 




1.4+0.2 


2-0 Q(l) 


1.238 


0.6±0.3: 


2.6+0.2 


2-0 Q(2)+4-2 S(4) 


1.242-1.242 




1.4+0.2 


2-0 Q(3) 


1.247 


0.9±0.4: 


2.8+0.2 


2-0 Q(4) 


1.254 




1.1+0.2 


4-2 S(3)+3-l S(0) 


1.261 




2.7+0.3 


2-0 Q(5) 


1.263 




2.1+0.3 


2-0 Q(6) 


1.274 




0.7+0.2 


4-2 S(2) 


1.284 




0.6+0.2 


2-0 Q(7) 


1.287 




1.9+0.2 


2-0 Q(8) 


1.302 




0.7+0.2 


4-2 S(l) 


1.311 




0.9+0.2 


5-3 S(5) 


1.312 




0.7+0.2 


3-1 Q(l) 


1.314 




1.2+0.2 


3-1 Q(2)+2-0 Q(9) 


1.318-1.319 




1.6+0.2 


3-1 Q(3) 


1.324 




1.3+0.2 


3-1 Q(4) 


1.333 




0.5+0.2: 


2-0 0(3) 


1.335 




2.7+0.2 


3-1 Q(5)+4-2 S(0) 


1.342-1.342 




1.6+0.2 


5-3 S(3) 


1.347 




0.7+0.2 


3-1 0(3)+4-2 Q(4) 


1.418 




1.3+0.4 


4-1 Q(5)+2-0 0(5) 


1.430-1.432 




2.1+0.4 



Position-velocity (P- V) diagrams are presented in Figure [6] 
(CGS4 Slits 1 and 2), and in Fig. [7] (CGS4 Slits 3 and 4). In 
Fig. [8] the emission line profiles are shown. It is worth to note 
that our spectra often cannot resolve the single structures inside 
the wiggling jet, due to the limited spatial resolution. As a result, 
radial velocities are often an average over more than one sub- 
structure visible in the Subaru image. 

Almost all the knots also have more than one velocity com- 
ponent in their spectral line profile, depending on the position 
of the slit with respect to the knot. However, such components 
are never completely resolved (with the exception of knot 2 
along the 'jet', see Fig. [8j, but they rather show up like 'bumps' 
along the smooth profile. When single-peaked profiles are ob- 
served, the lines are never symmetric, but evidence for extended 
line-wing emission (opposite to the blue- or red-shifted peak) is 
always detected. This may be indicative of a bow-shock mor- 
phology of the knots, producing both the wings and the differ- 
ent velocity components (see e. g. Davis et al. 120011 Schultz et 
al. 120051) . as also indicated by the sub-millimetric observations 
of Su et al. ( 120071 1. Remarkably, knot B has four velocity com- 
ponents (see Tab.|4j and Fig. [8}. The fourth component located 
at ~64kms~ 1 , and visible along the profile as a 'bump' on the 
red-shifted wing (see also knot X profile), cannot be explained 



in a bow-shock context. This component could originate from a 
different flow, but we have no further evidence in support of this 
hypothesis. 

The flux peaks of the knots range from -42 to -14kms _1 in 
the blue lobe, and from -8 to 47 km s _1 in the red lobe (see also 
Tab.|4]i. Noticeably, knot X and A in the red lobe (Fig.|6l central 
panel) have a slightly negative peak velocity. This can be under- 
stood considering that close to the source the axis of the flow has 
an inclination of ~9° with respect to the plane of the sky, and the 
aperture angle of the precessing jet is ~37°. Consequently, even 
if the knots are located in the 'red' lobe, they could have a nega- 
tive v ra d (see also Su et al. 200Tb. 

In both lobes, the absolute peak radial velocities of knots 
close to the source (A- • -C) are smaller (0-30kms _1 ) than 
those located at larger distances (knot D and knots 1 • ■ -4, 40- 
50kms _1 ). On the other hand, the FWZI of the line profiles de- 
creases with distance (see Tab. |4]and Fig. |8). The line profiles 
of the knots close to the source are broader on average ( 1 10— 
180 kms -1 ) than those far from the source (70-90 kms _I ). This 
could confirm that the inclination of the flow axis (with respect 
to the plane of the sky) is different in the two regions (i. e. it 
changes from ~9°, close to the source, to ~45° further out, at 
knots !• • 4 and knot D). On smaller scale, radial velocities ap- 
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Table 2. Observed lines in the IRAS 20126+4104 outflow: knots A, B, C, X, jet-1, jet-4. (Continued) 



1 



Term 








F ± A F(10" 


l3 ergcm 2 s 1 


) 




H L' 

2 — Mfjj 




A 




Q 




■ 

Knot 1-jet 


Knot 4-jet 


2-0 0(6) 


1.487 






1.0±0.4: 








3-1 Q(13) 


1.502 






0.6±0.3: 








3-1 0(5) 


1.522 






1.3±0.4 








2-0 0(7) 


1.545 






0.7±0.3: 








4-2 0(5) 


1.622 






0.9±0.3 








3-1 0(7) 


1.645 






0.7±0.2 








1-0 S(ll) 


1.650 






1.0±0.3 








1-0 S(10) 


1.666 






1.1 ±0.3 








1-0 S(9) 


1.688 


1.9±0.4 


2.1±0.3 


7.4±0.3 




0.9±0.3 


0.7±0.3: 


1-0 S(8) 


1.715 


1.4±0.4 


1.2±0.2 


5.2±0.3 




0.7±0.3: 


0.7±0.3: 


1-0 S(7) 


1.748 


12.1±0.4 


7.5±0.2 


33.7±0.3 


1.9±0.2 


3.8±0.3 


2.5±0.3 


1-0 S(6) 


1.788 


7.9±0.4 


5.0±0.2 


22.8±0.3 


1.4±0.2 


2.4±0.3 


1.5±0.4 


2-1 S(9) 


1.790 






1.0±0.5: 








2-1 S(8) 


1.818 






1.0±0.5: 








1-0 S(5) 


1.836 


54±5 


31±5 


125±5 








2-1 S(7) 


1.853 






8.0±5: 








1-0 S(4) 


1.891 


28±5 


16±5: 


66±5 








2-1 S(7) 


1.945 






12.0±5: 








1-0 S(3) 


1.958 


24±2 


51±5 


110±5 


13±4 






3-2 S(7) 


1.969 






1.7±0.8: 








2-1 S(4) 


2.004 




1.3±0.4 


8.3±0.8 




0.6±0.3: 




3-2 S(6) 


2.013 






1.2±0.5: 








1-0 S(2) 


2.034 


38.9±0.5 


21.1±0.4 


83.4±0.5 


4.2±0.4 


7.7±0.3 


5.7±0.4 


3-2 S(5) 


2.066 


0.9±0.3 


0.7±0.3: 


2.9±0.4 








2-1 S(3) 


2.073 


7.4±0.4 


4.8±0.3 


24.4±0.4 


1.4±0.4 


2.0±0.4 


1.1 ±0.3 


1-0 S(l) 


2.122 


121 ±0.4 


66.4±0.4 


246±0.4 


11.1±0.4 


22.6±0.4 


16.6±0.4 


2-1 S(2) 


2.154 


3.7±0.4 


2.2±0.4 


8.4±0.4 


0.5±0.2: 


1.0±0.3 


0.6±0.2 


3-2 S(3) 


2.201 


1.6±0.4 


0.8±0.3: 


4.3±0.4 


0.4±0.1: 


0.5±0.2: 




1-0 S(0) 


2.223 


30.1±0.4 


15.6±0.4 


59.5±0.4 


2.8±0.3 


5.7±0.3 


4.6±0.4 


2-1 S(l) 


2.248 


10.9±0.4 


5.5±0.4 


23.0±0.4 


1.6±0.4 


2.5±0.3 


2.1±0.4 


3-2 S(2) 


2.286 


0.9±0.4: 




1.5±0.4 








2-1 S(0) 


2.355 


2.7±0.5 


1.7±0.5 


5.1±0.5 








3-2 S(l) 


2.386 


1.0±0.5: 




4.0±0.5 








1-0 Q(l) 


2.407 


135±5 


66±5 


234±5 


11±3 


21±5 


16±5 


1-0 Q(2) 


2.413 


50±5 


25±5 


89±5 








1-0 Q(3) 


2.424 


123±5 


57±5 


218±5 


11±3 


20±5 


13±5: 


1-0 Q(4) 


2.437 


35±5 


18±5 


70±5 








1-0 Q(5) 


2.455 


39±5 


27±5 


142±5 








[Fe II] lines 


a 4 D lr _ -a 6 D 9l2 


1.257 




0.6±0.3° 




1.2±0.3 






a 4 D 7/2 - a 6 D 1/2 


1.321 








0.6±0.3° 






a 4 D 1/2 - a 4 F 9/2 


1.644 




1.7 ±0.2 




2.0±0.2 






a 4 D 1/2 - a 4 F ll2 +a 4 P 5 / 2 - a 4 D 1/2 


1.810-1.811 




1.2±0.4 




1.5±0.4 







Notes: 3o~ upper limits (in knot C) for the [C i] doublet (at 0.983-5 fim) and for the Bry line are 10" 15 and6xur 16 erg cm s , respectively. 



pear to oscillate. This is well visible in P-V diagrams of Slit 3, 
and marginally, for the red-shifted knots in Slit 1. 

Assuming two different values for the inclination, the spatial 
velocity of the knots ranges between 50 and 80 km s _1 . 



3.4. Physical parameters of the gas 

The wealth of molecular hydrogen lines detected in both MIR 
and NIR spectra allows us to perform a detailed study of the H2 
excitation in a high-mass jet, taking into account, for the first 
time, also the pure rotational lines. 

As a result, employing all the available H2 ratios in our ro- 
vibrational diagrams, we have derived column densities, extinc- 
tion, and temperature of the gas. Combining these parameters 
with the 2.12//m flux obtained from the narrow-band imaging, 
we have determined an accurate measurement of the H2 lumi- 
nosity (Lf/i) for each knot and for the entire flow (see Caratti o 
Garatti et al. 2006 for a detailed description of this procedure). 



3.4.1. Ro-vibrational diagrams from NIR lines 

As a first step of our analysis, we employed only the NIR lines. 
Line pairs originating from the same energy level should lie in 
the same position of the ro-vibrational diagram: by varying the 
extinction value (A v ) and increasing the goodness of the fit (max- 
imising the correlation coefficient) extinction and temperature 
can be evaluated simultaneously (see e. g. Giannini et al. 2004 
Davis et al. l2004l Caratti o Garatti et al. 2006). In order to reduce 
the uncertainties, only the transitions with a S/N>3 and not af- 
fected by blending with other lines were used. The result of this 
analysis is shown in Table [5] where the average excitation tem- 
perature, the A v , and the column densities of the warm gas com- 
ponent are listed for each knot (columns 2^1). The ro-vibrational 
diagrams are shown in Fig. [9] 

To compute the extinction in each knot, we have se- 
lected all the available pairs of lines. From the spectrum 
of knotC 41 lines have been used, giving 29 different ra- 
tios (i. e. the ratios are among the following groups of transi- 
tions: l-0S(i),l-0Q(i+2); 2-0S(i),2-0Q(i+2),2-0O(i+4),2-lS(i); 
3-lS(i),3-lQ(i+2),3-10(i+4); with i=0,l,2...). The best fit was 
obtained excluding the l-0S(i)/l-0Q(i+2) pairs, that deviated 
more than 3<x from the average value. The final value is 
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Fig. 4. 0.9-2.5 yizm low resolution spectrum of knot C in the IRAS 20126+4104 jet. An asterisk near the line identification marks 
the detections between 2 and 3 sigma. Telluric lines are indicated by the symbol "ffi". 



A v -1. 6+0.2, and has a very small error. For the other knots, the 
A v ranges from 6 to 10 mag, with errors between 1 and 3 mag 
(see Tab. [5] column 3). For these knots, a few pairs of lines were 
available and hence the errors are larger. 

With the exception of knot C, all the lines in our ro- 
vibrational diagrams are well fitted by a single straight line (see 
Fig. [9}, indicating a uniform temperature of the gas. The ob- 
served temperatures range from 2000 to 2500 K. 

On the other hand, the different H2 lines detected in knot 
C cannot be fitted by a single line. The ro-vibrational diagram 
exhibits a typical curvature indicating the presence of a strat- 
ification in the gas temperature. A more elaborate model of a 
mixture of gas at two different temperatures can describe these 
data (see e. g. Giannini et al. 2002 , Caratti o Garatti et al. 2006 
Gredel 2007 ). The population densities of the lines coming from 
levels with an excitation energy up to ~12000K are consistent 
with an excitation temperture of 2050 K. These lines indeed trace 
the warm component of the gas. Conversely, lines with a higher 
excitation energy (>12 000K) are thermalised at a higher tem- 
perture of -5200 K. These lines trace the hot component of the 
gas. A single fit through all the lines gives only a measure of the 
'averaged' temperature (~3300K) (see Fig. [9] top left panel). In 
our model the hot component is a fraction of the gas of about 
8%. 

Another important parameter that derives from ro- 
vibrational diagrams is the column density of the gas. The 
column densities of the warm H2 component are reported in 
Table [5] They are very similar to those measured in other 



high-mass protostellar jets (10 18 -10 19 cirT 2 ) (see e.g. Davis et 
al. 120041 Gredel 12006 ), and they are 1-2 orders of magnitude 
larger than the values observed in low-mass jets. 

Finally, from the analysis of Fig. |9j it is worth to note that 
the H2 gas along the flow is fully thermalised. The fact that 
we detect no lines with v>6 indicates that fluorescence mech- 
anisms do not play an important role in the excitation (see e. g. 
Black & van Dishoeck 1987). Excitation by non-thermal pro- 
cesses like fluorescence in presence of a UV field would cause 
the vibrational temperature to be different from the rotational 
temperatures, and thus lines belonging to different vibrational 
series should not be located on the same line (see e. g. Hora & 
Latter [19941 ) . As a consequence, we should observe strong de- 
viations from the smooth (or linear) Boltzmann distribution of 
our ro-vibrational diagrams (see also Gredel 2007). Apparently, 
this contrasts with results of Sec. 13.2.21 where the presence of a 
PDR was inferred. Such a PDR could produce enough far-UV 
radiation to excite the higher vibrational levels of the H2 and 
induce a fluorescent emission (see e. g. Burton 1992 ), that we 
do not detect. To evaluate a possible contribution of the PDR to 
the observed lines, we used the 'PDR toolboxfl The observed 
[Cn] 158//m diffuse line intensity (considering a LWS beam- 
width of 80") is ~2.5xl0 _4 ergs~ 1 cirT 2 sr _1 . Assuming that all 
the emission arises from the PDR and that the density is 10 5 - 



3 The 'PDR Toolbox' is available at http://dustem.astro.umd.edu and 
contains downloadable FIR line diagnostic information about PDRs. 
The tool has been created by L. Mundy, M. Wolfire, S. Lord, and M. 
Pound, and it is based on the new PDRmodels of Kaufman et al. (1999). 
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X Offset (arcsec) velocity (km/s) velocity (km/s) 

Fig. 6. P- V diagram for slits 1 and 2. The left panel shows the H2 (TNG) image with the slits superimposed, and the position of 
the knots. Y-offsets are displayed from the position of the source. The contour levels are 3, 30, 50, 100, 130, 210, 235, 300, 420x 
the standard deviation to the mean background (<x ~ 2 x 10~ 17 erg s _1 cirT 2 arcsec -2 ). The central and right panels display the 
uncalibrated H2 (2.12/mi) spectra from CGS4 (slits 1 and 2, respectively). The contour levels are 5, 20, 60, 125, 180, 300, 450x the 
standard deviation to the mean background. 
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Fig. 8. Line profiles of the knots in the 1-0 S(l) line of H2 at 
2.12 yizm. 



10 6 cm -2 (Cesaroni et al. 1999), we obtain a relatively faint FUV 
flux of Go between 10 and 100 (where Go is measured in units 
of 1.6xlO -3 erg s _1 cm 2 ). Under these circumstances the PDR 
surface temperture does not exceed 100 K (Kaufman et al. 1999 ) 



and the contribution to the overall H2 emission is negligible. For 
example, the diffuse line intensity of the 0-0 S(l) is between 10~ 6 
and 10~ 5 erg s _1 cirT 2 sr _1 , while from our observations we ob- 
tain ~6xl0~ 4 erg s~' cirT 2 sr _1 (dividing by the ISO-SWS FoV). 
In the NIR the emission from the shock (at 2.12/im) is four or- 
ders of magnitude larger. These results clearly indicate that the 
detected H2 mostly arises from shocks. 

3.4.2. The cold H 2 component from MIR lines 

As a second step in our analysis, we also used the ISO-SWS lines 
to study the cold component of the gas, that is usually traced by 
lines with excitation energy lower than 5000 K, namely the 0- 
0S lines (with J fi „ <5) between 6 and 28jum in the ISO-SWS 
spectrum. 

Since we do not know the spatial extent of the emitting re- 
gion, we have matched the column densities of the 1-0Q lines 
present in both NICS and SWS spectra, to intercalibrate the data 
in the ro-vibrational diagram. Indeed the measured ISO-SWS 
fluxes of the 1-0Q lines are very close to those measured in 
knotC. Therefore we have assumed that its column density is 
representative of the warm H2 component detected with ISO, as 
well. The result is shown in Fig. [10] SWS lines coming from the 
v=0 and v= 1 levels are represented as open and filled pentagons, 
respectively. The data are overplotted on the ro-vibrational dia- 
gram of knot C (shown in Fig. |9j top left panel). Although there 
is some scatter among the v=l data points (two 1-0Q and l-OO 
lines have slightly larger and smaller column densities, respec- 
tively), the overall fit is quite satisfactory. In particular, the 0- 
S lines reveal the presence of a cold gas component at about 
520 K, that has a larger column density than the other compo- 
nents (N H , ~9.7xl0 21 cm~ 2 , see also Fig.fTOl. 

In order to reproduce the observed ro-vibrational diagram, 
we calculated a theoretical H2 spectrum for a mixture of three 
H2 layers in LTE condition at a temperature of 520, 2050, and 
5200 K. The adopted LTE code (see also Caratti o Garatti et 
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X Offset (arcsec) velocity (km/s) X Offset (arcsec) velocity (km/s) 

Fig. 7. Left P-V diagram for slit 3. On the left panel the H2 (TNG) image indicates the position of the slit and knots. Y offsets are 
displayed from the position of the star in the centre of the image. The contour levels are 3, 6, 12, 20, 30, 50x cr. On the right the 
uncalibrated H2 (2.12/im) spectrum from CGS4 is reported. The contour levels are 3, 5, 7, 10, 13x cr. Right P-V diagram for slit 
4. Y offsets are displayed from the position of knot D. The image contour levels are 3, 6, 30, 50, 100, 130, 210, 235, 300, 420x cr. 
The levels of the spectra are 3, 5, 7, 10, 13x cr. 



al. [2006) computes the line intensities involving levels with 
<v< 14 and < / < 29 (E,, y < 50 00 K). Ro- vibrational 
energies are taken from Dabrowsy (1984) and the Einstein co- 
efficients from Wolniewicz et al. (1998 ). We have assumed an 
ortho/para ratio equal to three (as also evinced from our ro- 
vibrational diagrams). The resulting model is plotted as a dashed 
line on the points of Fig.[l0] and properly fits our observations. 
In this case the warm and hot components are only a small frac- 
tion of the total H2 gas, i. e. less than 1%. 

3.4.3. H 2 luminosity 

In the last step of our analysis we inferred . Once the physical 
parameters (A v , T) for each knot were derived, we dereddened 
the 2.12/im flux obtained from the imaging (see Tab. [5] column 
5), adopting the Rieke & Lebofsky (1985 ) reddending law, and 
computed the line ratios with the other H2 lines by applying our 
radiative LTE code at T = T avg . From these ratios, the abso- 
lute intensities of individual lines were computed and the H2 
luminosity of each knot was derived (see Tab. [5] column 6). The 
sum of these values gives a total of the flow of 4.6+0.3 Lq. 
An average temperature reproduces the contribution of both the 
warm and hot gas components to the radiated energy in H2 well 
(see Caratti o Garatti et al. 2006 ). We are not taking into account 
the luminosity of the 'cold' component, that, however, due to the 
extremely high column density of the gas, here can be consider- 
able. In order to evaluate this, we have assumed that the mea- 
sured SWS flux from the 0-0 S lines is the total flux emitted by 
the shocked cold H2 gas of the flow, and that the contribution of 
non-thermal emission is negligible, as seen in Sec l3.4.ll We then 
computed the luminosity of an LTE gas at T=520K, using the 



0-0 S(l) observed dereddened intensity to derive the absolute in- 
tensities of the remaining lines. The luminosity of the 'cold' gas 
is 3.6+0.6 Lq. Added to the previous estimate this gives a total 
Lh 2 for the entire flow of 8.2+0.7 Lq. 

3.5. Mass flux, mass, and energy of the H 2 jet 

The physical and kinematical parameters inferred in the previous 
sections allow us to evaluate important kinematical and dynam- 
ical properties of the flow. 

Therefore, we have estimated the mass flux rate of the flow 
(M ol „) from the H2 in order to compare it with the mass outflow 
rate, and the mass accretion rate of the protostar previously ob- 
tained from CO observations and models (see e. g. Shepherd et 
al. [20001 Cesaroni et al. [20051 Lebron et al. [2006l >. 

M 0UI (H2) can be written as Mj. = 2fim H NH^Av,/l h where yu 
is the average atomic weight, thh the proton mass, Ng 2 the H2 
column density, A the area of the H2 knot, v, the tangential ve- 
locity, and I, the projected length of the knot (see, e. g., Nisini et 
ai. |20"0"5l Podio et al. 120061 Antoniucci et al. |2007l >. 

Using the 2.12 fim column density of Tab. [5] the radial ve- 
locities of Tab. [4] and assuming an average inclination for all 
the knots of 20° with respect to the plane of the sky, we esti- 
mate mass fluxes between 10~ 6 and 10~ 8 MQyr~' (see Tab. [5] 
column 7). In fact, these values represent lower limits. From the 
ro-vibrational analysis the cold gas column density is two orders 
of magnitude larger than the warm component (e. g. 0-0 S lines 
have a column density Nh 2 ~9.7x10 21 cm -2 ). Using this value, 
and assuming the same velocity and extension of the emitting 
area of the 2. 12/vm line, we obtain M ou ,(H 2 )~7.5x 10~ 4 M Q yr _1 . 
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Fig. 9. Rotational diagram of knots in the IRAS 20126+4104 jet obtained from the low resolution NIR spectroscopy. Different 
symbols indicate lines coming from different vibrational levels, as coded in the upper right corner of the boxes. The infeiTed 
temperature and extinction is indicated in the upper right corner of the box. 



It is worth noting, however, that this value should only be 
considered as an estimate of the mass flux, for two reasons. 
Firstly, we have set I, and v, equal for both H2 components. 
Probably, such an estimate is a lower limit, since l t of the cold gas 
can be larger than the I, measured at 2.12/im. More important, 
the adopted mass flux formula lays on the assumption that l t /v t 
is the cooling time (f c ) of the shock. In this case, the kinemati- 
cally determined t c of knot C is 1.9xl0 1() s. Also other methods 
have been used in the literature to derive t c : Davis et al. (2000) 
adopt the approximation t c ~ 3 X 10 8 n 6 1 rr 2 3 (in seconds), 
where «e is the density in units 10 6 cirT 3 and T3 is the tem- 
perature in units of 1000 K (Smith & B rand 1 19901 . In this way, 
adopting n# 2 = 10 5 cirT 3 (Cesaroni et al. 119991 ) and two different 
temperatures, for knot C, we would obtain f c (warm)= 6.1x10 s s 
(with r=2000K), and * c (cold)= 1.5xl0 10 s (with r=500K). The 
M out of the cold component would then be remarkably similar 
to the previous estimate. Moreover, it is very close to the mass 
outflow rates derived by Shepherd et al. (f2000) and Lebron et 
al. d2006b from the CO mm analysis (8.1xlO" 4 M Q yr" 1 , and 
34xl0~ 3 M© yr -1 , respectively). 

In addition, we can derive a further estimate of the mass flux 
from the luminosity of the [0 1] line (at 63 /mi) using the approx- 



imate formula M out ~10 4 x(L(63/mi)/L )MQyr 1 (see e.g. 
Hollenbach & McKee |T989l Liseau et al. flWTl Cabrit |2002l ). As 
a result, we get M out (OI) = 2xlO~ 4 M0yr -1 . Such an estimate 
is, however, an upper limit, since part or all of the [0 1] emission 
could originate from the PDR. 

To derive the mass of each knot, we assume Mk - 
2/umHniHNHnA. For knot D, where Nh 2 was n °t measured, we 
have assumed an average value of 4.4xl0 18 crrT 2 from Tab. [4] 
The total mass (S^Mt) of the warm gas is ~10~ 3 Mq. Assuming 
an average value for the inclination of the flow of 20°, the 
dynamical timescale of the flow (r^) is ~1.3xl0 4 yr, the to- 
tal momentum PCE^MkVk) is 0.06Mokms _1 , the kinetic en- 
ergy Ejt(iSjfcMjfeV?) is 5xl0 43 ergs, and the momentum flux P 
(P/ T d) is 5xl0~ 6 Mq yr~' km s~' . We note that these values must 
be considered as lower limits, since the mass of observable 
shocked warm H2 is only a small fraction of the total shocked 
molecular hydrogen mass, as seen previously. If we compute 
these quantities taking into account the column density of the 
cold component of the gas, we obtain the following crude esti- 
mates: M ~O.6M0kms _1 , P -SOMokms" 1 , E k ~4xl0 46 ergs, 
p ~4xlO- 3 M yr- 1 kms- 1 . 
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Table 5. Physical parameters of the H2 knots in IRAS 20126+4104 derived from the low resolution spectroscopy and imaging. 



Knot 


Average Temperature 


A V (H 2 ) 


Nh 2 


Flux (2. 12 /an) 




M out (H 2 f 




(K) 




10 l8 cnT 2 


(10~ 14 ergcnT 2 s~') 




(10~ 8 Mrr, vr" 1 ) 
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10+3 
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29.9+0.2 


0.70±0.13 
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£-1-1 c 
0± 1 


3.0 


1 n Q-t-f» 1 


U.Zl ±U.UJ 
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3300+100 


7.6+0.2 


12.5 


71.5+0.3 


3.12+0.24 


192 
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6.4+0.1 


0.1 2 d 
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2500+100 


6±l f 


1.9 


2.4+0. l f 


0.06+0.01 
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jet - knot 1 


2350+90 


6+3 


1.7 


7.5+0.2 


0.2+0.1 


8 


jet - knot 2 








1.0+0.4 


0.02 s 




jet - knot 3 








2.7+0.5 


0.04'' 




jet - knot 4 


2020+70 


6+3 


1.2 


10.0+0.2 


0.2+0.1 


8 



Notes: " Derived from the average temperature. The contribution from the cold component of the gas is not taken into account. The L Hl from the 
cold component, as measured from the 0-0 S lines of the ISO spectrum, is L Hl =3.6+0.4 Lq. 

* Derived from the average temperature. The total N m column density is measured by comparing the intrinsic (per unit mass) and observed flux 
ofthe 1-0 S(l) (2.12 /mi) line. 

c A v measured from the [Fen] lines (1.64 and 1.26 /mi) is 11+7 mag. 
d Computed assuming r„=2000K and A v =6 mag. 
e A v measured from the [Fe 11] lines (1.64 and 1.26 /mi) is 7+4 mag. 
f Flux([Fen]) = 3±lxl0~ I5 ergcrrr 2 s~' 
s Computed assuming T„ =2350 K and A v =6 mag. 
Computed assuming T„=2020 K and A v =6 mag. 



Finally, we can give a very rough estimate of the density 
for the cold component, as well. If we assumed that the size for 
knot C along the line of sight was similar to the observed length 
(~10 17 cm), we would obtain ng, ~10 5 cm~ 3 , as in Cesaroni et 
al. (fT999b . 



4. Discussion 

4.1. A precessing jet model 

As mentioned in Sect. I3.ll the inner region of the H2 jet clearly 
shows a wiggling morphology, possibly indicating a small-scale 
precession, different from the long-period precession derived 
from the positions of the outer knots (Shepherd et al. 2000, 
Cesaroni et al. |2005l >. 

In order to investigate this new morphological feature, we 
used a simple jet model to fit the observed position of the inner 
knots (see e. g. Eisloffel et al. 1996, Texeira et al. 2007): 



a 



a 0\ + l cos{ijj) -sin(tff) \ I (plsin(2nl/A +Xo) . ([ 
6 I \ sin(tff) cos(i//) j \ I * 



where ao, So are the coordinates of the source, if/ the rotation an- 
gle, <p the precession amplitude (in radiants), / the distance from 
the source position, A the precession length scale, xo the initial 
phase at the source. The resulting model, shown in Fig.Qj] has 
been obtained for values of ^= -60!9, 0=0.133 (~7!6), A -1079, 
;to=8 1°.6. The inclination of the flow with respect to the plane 
of the sky was not considered in the fit, since its value (close to 
the source) is relatively small (~9°). For a constant jet velocity 
of 80kms _1 , at a distance (D) of 1.7kpc, we obtain a preces- 
sion period of ~1100yr, that appears short in comparison with 
the value of -64000 yr derived in Shepherd et al. ( 120001 ) and 
Cesaroni et al. (2005). We have then tried to reproduce the com- 
plex motion of the flow, adding a second precession mode to our 
model, that should take into account the longer period traced by 
the outer knots. We were, however, not able to fairly match the 
pattern of the entire jet. 



This could indicate that the origin of the two observed 'pre- 
cessions' is different, or that the precession period changed over 
time. The long-period precession was interpreted by Shepherd et 
al. (|2p00) and Cesaroni et al. (2005) as caused by the interaction 
between the disc of IRAS 20126+4104 and a stellar companion 
of a few solar masses. However, it is unlikely that the small- 
period precession we observe is due to the orbital motion of the 
jet source around its companion. Since the orbital radius is given 
by the following equation (see Anglada et al. 120071 1: 



Atan((f))D 
2n 



(2) 



the resulting value is r ~400AU. The radius of 
IRAS 20126+4104 disc was es timated around 800 AU on 
the contrary (Cesaroni et al. 1997), meaning that the companion 
would intersect the disc. Alternatively, the precession could 
have changed with time. The ejection of a third companion in 
a hierarchical triple system, for example, might have led to 
the formation of a tighter binary system and to a shortening of 
the orbital period. However, this hypothesis does not explain 
the different position angles on the sky measured in the two 
precession modes (-60°9 and -37°, respectively). 

It seems more likely, then, that the complex axis wandering 
is induced by tidal interactions of multiple stellar companions. 
Indeed IRAS 20126+4104 is not a single compact source, but 
it is composed by a small cluster of YSOs (De Buizer 120071 1, 
therefore the dynamics of such a system can be extremely com- 
plicated. It is indeed beyond of the scope of this paper to obtain 
a rigorous modelling for the dynamic of such a system. 

4.2. H 2 jet vs CO outflow 

Observations of protostellar jets in low-mass YSOs support the 
idea that their molecular outflows are driven and accelerated by 
the highly-collimated bipolar jets. Outflows from intermediate- 
and high-mass sources, on the other hand, are often poorly colli- 
mated and only in a few cases there is clear evidence of a driving 
jet. This has often led to the suggestion that these outflows could 
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ov=0 T=520±70 K 
ov=l T=2340±70 K 
. v=2 T = 3590±250 K 
£,v=3 T=5500±500 K 
v=4,5 T=5100±700 K 

- T, + T 2 + T 3 = 

(520 + 2050 + 5200) K 

A„= 7.6±0.2 mas 




Table 6. Comparison between H2 jet and CO outflow physical 
properties of the IRAS 20126+4104 flow. 
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Fig. 10. Rotational diagram of knot C of the IRAS 20126+4104 
jet obtained from the low resolution spectroscopy. Different 
symbols indicate lines coming from different vibrational levels, 
as coded in the upper right corner of the box (v=0 and v=l ISO 
lines are indicated by open and filled pentagons, respectively). 
The dashed line represents the theoretical population distribution 
for a combination of thermalised gas at three different tempera- 
tures (520, 2050, and 5200 K, see text). The derived temperature 
from the different vibrational levels is indicated in the upper right 
corner of the box. The extinction is reported as well. 




Fig. 11. Precession model plotted over the H2-Bry image from 
Subaru. The model is for a precession angle of 7°6 and a preces- 
sion scale of 10f.'9. A cross indicates the position of the source. 



be driven by a different mechanism, as, e. g., a wide-angle radial 
wind (see e. g. Arce et al. |20071 >. 

Apparently, IRAS 20126+4104 represents a particular ex- 
ample of poorly collimated outflow powered by a protostellar 
jet. The poor collimation of this outflow has been explained by 
means of the severe precession of the jet. It is thus interesting to 
investigate if the properties of the IRAS 20126+4104 H 2 jet meet 
the physical parameters of the CO outflow, i. e. if the jet is pow- 
erful enough to accelerate the surrounding medium and drive the 



Parameter 



IP 



CO" 



CO" 



Mass (Mo) 
v red (kms -1 ) 
\ blm (kms- 1 ) 
T rf (10 4 yr) 
M(10^M o yr) 
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/ , (lQ~ 3 M yr'kms~' 
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53 
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(4, 55) 
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0.75 


0.81 


3.4 


4 


5.1 
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50 
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1490 


4 


6 
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Notes: 

" From Shepherd et a l. (|2000l 
b From Pebron et al. ( I20061 . 



outflow. This is indeed a fundamental question, because the YSO 
accretion rate is often estimated from the mass flux rate of the 
outflow. Comparing CO and H2 parameters, however, we have 
to keep in mind that the CO lines only give a time-integrated 
response over the lifetime of the jet, whereas the H2 emission 
provides an 'instantaneous' measure of these quantities. 

Parge-scale outflow properties from CO mm observations 
were extensively studied by several authors (see e. g. Shepherd 
et al. 120001 Pebron et al. 120061 1. Since their computations are 
based on different assumptions, however the results of Pebron et 
al. (2006 ) are at least one order of magnitude larger than those 
of Shepherd et al. (120001 ). 

In Table [6] we compare our H2 flow parameters (obtained 
including the cold component) with the parameters of the CO 
outflow from Shepherd et al. d20001 > and Pebron et al. (2006). 
From this table, we can infer that the jet is, at least partially, 
driving the outflow. The radial velocities of H2 are well consis- 
tent with the CO velocities, including the high velocity compo- 
nent up to 60km s _1 detected by Pebron et al. (2006). On the 
other hand, the total mass of the Fp bullets is just a fraction of 
the total CO mass. As a consequence, the momentum (P) of the 
jet is at least one order of magnitude smaller than that of the 
outflow. However, our estimates of the mass flux, kinetic en- 
ergy, and momentum flux are, indeed, almost coincident with 
Shepherd et al. (2000), supporting the hypothesis that the Fp jet 
is fully driving the flow. Contrarily, the values from Pebron et 
al. ([2006 ) are from one to two orders of magnitude larger than 
ours (see Tab. [6]). This would imply that the Fp jet is only par- 
tially driving the outflow, and that there is an undetected compo- 
nent along the flow, likely a neutral jet, that is supplying further 
momentum flux to the outflow. However, it is worth to note that 
the mass flux obtained by Pebron et al. (2006) is slightly larger 
than the mass accretion rate of IRAS 20126+4104 inferred by 
Cesaroni et al. ( 120051 ) (2xl0~ 3 Mq yr" 1 ), and this is quite un- 
likely. This discrepancy could be then explained by the presence 
of (undetected) multiple jets (winds), driven by other YSOs near 
the source, that would substantially contribute to power the CO 
outflow. 



4.3. Shock conditions along the jet 

As the previous results point out, the cold H2 component of the 
jet plays a major role in the kinematics and dynamics of this out- 
flow. It is therefore desirable to discuss the origin of this compo- 
nent in more detail. 
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A prevalence of C-type shocks along the flow could explain 
our findings. Indeed C-type shocks produce a high column den- 
sity in the 0-0 lines (see e. g. Flower et al. 120031 McCoey et 
al. [2004), that are very sensitive to these shocks. The column 
densities in our 0-0 lines could not be explained by a J-type 
shock (see e. g. Le Bourlot et al 120021 1. Moreover, except for the 
[Fen] and the [Oi] emission close to the source, the emission 
observed along the flow comes only from the H2. Only the inner 
region near the source position exhibits hints of ionic emission. 
There is no evidence of the ionised jet detected in the radio by 
Hofner et al. (2007). probably confined to a region very close 
to the source and highly extincted. On the other hand, the com- 
plete absence of any ionic emission in the other knots of the flow 
should not be attributable to a high value of extinction, since the 
measured A v values along the jet are relatively small. Such a be- 
haviour has already been observed in several low-mass protostel- 
lar jets. They usually show only H2 emission, or, at most, a hint 
of ionic emission close the exciting source (see e. g. Giannini et 
al. 120041 Caratti o Garatti et al. |20061 >. It is also well observed in 
the high-mass IRAS 18151-1208 jet spectra (Davis et al. l2004b . 
They are interpreted as C-type shocks, or as J-shocks with mag- 
netic precursors, that is J-type shock waves evolving into "con- 
tinuous" C-type shock waves (see e. g. Giannini et al. [2004). 
Also our kinematical observations seem to indicate the presence 
of C-type shocks. H2 velocities in C-type shocks can exceed 
50-60 km s" 1 up to 80km s _1 (see e.g. Le Bourlot et al. I20021) . 
Accordingly, the FWZI of H2 bow-shocks can be up to twice this 
value, depending on the inclination of the jet axis with respect 
to the plane of the sky (see e. g. Davis et al. 120041 1. that is the 
broadest lines are expected for bow-shock located in the plane 
of the sky. Actually, if we consider the different inclinations for 
the inner and the outer knots, our observations seem to indicate a 
high shock velocity in both regions, whereas the larger value of 
the FWZI is observed in the inner knots (closer to the sky plane), 
and a less broad shape in the outer knots. From these values we 
can roughly estimate a shock velocity between 40 and 80 km s , 
where the fastest speed would be attained in the inner region. 
Such high shock velocities (50-80 km s~') however imply low 
pre-shock densities of 1 3 - 1 4 cm~ 3 (Le Bourlot et al. 2002} in C 
shocks, assuming a magnetic field Z?(/iG)=b[nH(cirT 3 )]~ - 5 , with 
b=l. Higher densities lead to a lower maximum shock speed 
(vrfiu), that is the maximum shock velocity which can be at- 
tained, prior to the collisional dissociation of H2. On the con- 
trary, both CO and H2 observations seem to indicate a higher pre- 
shock density (10 5 -10 6 cirT 3 ) in the medium (see e.g. Cesaroni 
et al. |1999) . In this case v^iss would decrease to 30-50km s _1 (Le 
Bourlot et al. 2002). Thus a 'standard' C-type model cannot en- 
tirely explain our findings close to the source. A larger transver- 
sal magnetic field (b>l) or a different shock model (as a J-shock 
with a magnetic precursor) are probably needed. This last hy- 
pothesis would also explain the presence of the ionic emission 
in the knots close to the source. 



4.4. Lh 2 , accretion and ejection rates 

Several authors have argued that a tight correlation between the 
evolutionary properties of YSOs and their outflows would ex- 
ist. This is particularly true for those objects forming by disc 
accretion, since accretion and ejection should be regulated by 
the same mechanism. Recently, it has been shown for low-mass 
YSOs that the total H2 luminosity of the jets is proportional to 
their accretion rates (see e. g., Smith 2002 , Froebrich et al. 2003 ). 
In particular, an empirical relationship between and Lb i (L^ 
oc L° h ^, Caratti o Garatti et al. 2006 ) was derived for a large sam- 
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Fig. 12. L H2 vs L bo i including the IRAS 20126 + 4104 jet (also 
considering the luminosity of the cold component, the new data- 
point is located in the upper-right corner of the diagram). Values 
for IRAS 18151-1208 and IRAS 1 1 101-5829 jets have been in- 
cluded as lower limits, as well (see discussion in the text). A 
dashed line indicates the previous fit from Caratti o Garatti et 
al. d2006l 



pie of low-mass protostellar jets. This correlation holds for very 
young YSOs (Class 0, and some Class I), where the bolometric 
luminosity is mostly coincident with the accretion luminosity of 
the object, and implies that M acc increases with the luminosity 
(i. e. mass) of the protostar. In Fig. [12] we report the results ob- 
tained from the sample of Caratti o Garatti et al. ((2(306'). The 
diagram compares the measured outflow H2 luminosity versus 
the bolometric source luminosity, both on a logarithmic scale. 
We have included the new data from IRAS 20126+4104 (posi- 
tioned in the upper right corner). The dashed line indicates the 
best fit previously obtained, that fits perfectly to our object, sug- 
gesting that the relationship applies also to more massive jets 
in their earliest stage of formation. Accordingly to the dynami- 
cal timescales of the outflow and the jet (few times 10 4 years), 
IRAS 20126+4104 has not yet reached the main sequence (MS), 
and has not yet developed any hypercompact HII (HCHII) re- 
gion, that may affect the collimation of the jet/outflow system 
(Beuther & Shepherd 2005). Most importantly, the bolomet- 
ric luminosity of the source mainly comes from accretion (see 
Cesaroni et al. 119991 Cesaroni et al. 2005). In these works, 
the authors obtain for this source an accretion luminosity of 
1.2x10 4 Lq for a mass accretion rate of ~l-2xl0~ 3 Mq yr _1 . 
Moreover, the fact that we obtain an ejection rate of 30-40% of 
the accretion rate from the H2 strengthens the reliability of our 
findings, since for high-mass sources we would expect a ratio 
of the mass ejection/accretion rate larger than the value (~0.1) 
usually derived for the low-mass objects (see e. g. Cabrit 2007 ). 

Finally, we compare our results with the L# 2 of two high- 
mass jets previously investigated by means of NIR spectroscopy, 
i.e. IRAS 18151-1208 (L Hl =QJ L Q , D avis e t al. 120041 ) and 
IRAS 11101-5829 (L Hl >2L Q , Gredel [2006t . These sources 
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have almost the same Lb i as IRAS 20126+4104. The Lh 2 esti- 
mates are quite close to the H2 luminosity of IRAS 20126+4104 
obtained from the warm H2 component (i. e. from our NIR 
analysis, Lh 2 =4.6+0.3 L0). For comparison, the values of the 
IRAS 18151-1208 and IRAS 1 1 101-5829 jets have also been in- 
cluded as lower limits in Fig. Q~2] It is also worth to note that 
the slightly lower value found by Davis et al. (2004) could be 
caused by the high extinction values (10-30 mag) observed to- 
wards this flow. Indeed a MIR investigation could reveal if the 
cold H2 component plays a major role in the cooling of those 
jets, as well. 

4.5. Comparing high- and low-mass jets 

The three high-mass protostellar jets spectroscopically inves- 
tigated up to date share similar characteristics. As the low- 
mass jets, they are collimated and powerful enough to drive 
their outflows. The jet is formed close to the source, as we 
observe e.g. in IRAS 20126+4104 knot X (located -1700 AU 
from the source), or in IRAS 1 1 101-5829 knot HH136J2 (lo- 
cated ~3000AU from the source), indicating that the collima- 
tion of the jet occurs close to the source. They show molec- 
ular and ionic shocked emission along the flow, and no evi- 
dence of fluorescent excitation is detected. The observed jet ve- 
locities are similar to those of the CO outflow. Where the ex- 
tinction is low, some jets (eg. IRAS 11101-5829, IRAS 18162- 
2048) show HH objects, as well. Some kinematical and dynam- 
ical quantities as mass flux, momentum flux, luminosity, and ki- 
netic energy, are larger however, than in low-mass jets, because 
the powering YSO is more massive. Moreover, a precession- 
wiggling morphology is observed in most of the massive col- 
limated H2 jets, often with large precessing angles (up to 
-4 0° in I RAS 20126+4104) (e. g. IRAS 1 6547-4 247, Brooks et 
al. 120031 IR AS 18 151-1208, Davis et al. 20041, IRAS 1 1101- 
5829, Gredel[2006] M17 disc silhouette, Niirnberger et al . 120071 
IRAS 07427-2400, I RAS 2 0293+3952, and IRAS 23033+5951, 
Nanda Kumar et al. 120021 and IRAS 20126+4104). Precession 
is also observed in low-mass jets, but it is not as frequent and 
as pronounced (usually the precession angles are less than 10°). 
All this indicates that the dynamical interactions among mas- 
sive stars are stronger and more frequent than in low-mass star 
forming regions. Such complex dynamics could also explain the 
confused H2 morphology of some jets and the lack of colli- 
mation in some massive outflows. Finally, the lack of jets and 
the poor outflow collimation observed in several massive young 
sources fits well the evolutionary outflow scenario proposed by 
Beuther & Shepherd d2005l l. In this context, those sources that 
have jets detected toward them are very young (well before the 
MS turn-on), while those without detectable jets near the pro- 
tostar have ultracompact HII (UCHII) regions. If the disappear- 
ance of a collimated jet in early B protostars is due to the pres- 
ence of enhanced ionising radiation from an accreting early main 
sequence star, then all early B stars may be formed via accretion. 
In this sense, they are scaled-up versions of low-mass protostars 
at early phases (as for IRAS 20126+4104, IRAS 18151-1208, 
IRAS 11101-5829, etc.). As YSOs evolve, developing UCHII 
regions that destroy the disc, the jets finally disappear, and the 
outflows look more like poorly collimated wind-blow bubbles 
(see also Arce et al |20071 >. 

In conclusion, the three high-mass protostellar jets spec- 
troscopically studied appear as scaled-up versions of the low- 
mass ones. Furthermore, a morphological analysis of the few 
intermediate- high-mass H2 jets known up to now partially sup- 
ports an evolutionary jet/outflow scenario. A larger sample of 



intermediate-/high-mass jets is, however, needed, and it would 
be premature jumping to the conclusion that the disc accretion- 
ejection paradigm can be extended to the intermediate- and high- 
mass protostars. 



5. Conclusions 

The IRAS 20126+4104 H2 jet has been extensively investigated 
through near-IR H2 and [Fen] narrow-band imaging, H2 high 
resolution spectroscopy, along with low resolution spectroscopy 
(0. 9-200 fim) throughout the infrared wavelength range. The 
kinematical, dynamical, and physical conditions of the H2 gas 
along the flow have been probed. The main results of this work 
are the following: 

- A high angular resolution H2 continuum subtracted image 
from Subaru reveals a small-scale precession of the jet close 
to the source, with an angle of (~7°6) and a period of 
~1100yr. This is about a factor 50 shorter than the preces- 
sion period deduced from large-scale H2 images. 

- H2 and [Fe 11] narrow-band images show the appearance of 
a new knot, labelled X, roughly 1" from the source posi- 
tion. No further ionic emission is detected along the flow in 
the narrow-band imaging, indicating that the jet is mainly 
molecular. 

- Low resolution spectra are rich in H2 emission, and no ionic 
emission is detected along the flow, with the exception of 
faint emission of [Fe 11] close to the source position (in knots 
X and B). Faint [Oi] and [Cn] emissions are observed in 
the ISO-LWS spectrum. They could arise from an embedded 
PDR region around the source. 

- The peak radial velocities of the knots range from -42 to - 
14kms _I in the blue lobe and from -8 to 47kms~' in the 
red lobe. Their line profiles, very broad and often with two 
or three velocity components, seem to indicate a bow-shock 
structure. In both lobes, the absolute peak radial velocities of 
knots close to the source (A- • C) are smaller (0-30kms _1 ) 
than those located at larger distances (i. e. knot D, and knots 
1 to 4, 40-50kms~ 1 ). This possibly confirms a change of 
the flow inclination angle (with respect to the sky) from ~9° 
(close to the source), to -45° (in the outer knots). Assuming 
these inclination values, the spatial velocity of the knots is 
between 50 and 80 km s~ 1 . 

- The ro-vibrational diagrams indicate H2 excitation tempera- 
tures between 2000 and 2500 K. Stratification of temperature 
is detected only in knot C, which can be modelled combining 
a warm (r ex =2050K) and a hot (~5200K) H2 component. 
Additionally, ISO-SWS spectrum reveals the presence of a 
cold component (520 K) with a high column density. 

- Furthermore, our analysis seems to indicate that the H2 is 
mostly excited in C-type shocks, and no evidence of fluores- 
cent excitation has been observed. 

- The estimated Lh 2 is 8.1 +0.7 L0, where the cold compo- 
nent contributes about 50% to the whole radiative cooling. 
The large H 2 luminosity suggests that IRAS 20126+4104 
has a significantly increased accretion rate compared to the 
low-mass YSOs. This is also supported by the measured 
mass flux rates from H 2 lines (M o „ f (H 2 )~7.5xlO~ 4 M0 yr _1 ) 
well matching the previous CO estimates. Our analysis also 
points out that the cold H2 component plays a major role in 
the kinematics and dynamics of this flow. 

- Comparing the H2 and outflow parameters strongly indicates 
that the jet is driving, at least partially, the outflow. 
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- By comparing the measured luminosity of the H2 jet with 
the source bolometric luminosity (assumed representative of 
the accretion luminosity), we show that IRAS 20126+4104 
fits well the correlation between these two quantities al- 
ready found for low-mass protostellar jets (Caratti o Garatti 
etal.|20p. 

- Considering our results and literature data of a few 
intermediate- and high-mass protostellar jets, we conclude 
that these few jets appear to be scaled-up versions of their 
low-mass protostellar counterparts. 
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